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1.  INTRODUCTION 


Acoustic  backscatter  from  the  seafloor  has  been  a  topic  of  considerable 
interest  over  the  past  several  decades.  The  importance  of  seabed  acoustics  has 
recently  grown  due  to  current  emphasis  on  mine  warfare  in  the  littoral 
environment.  Since  there  is  at  present  no  general  agreement  on  the  physics  of 
the  scattering  process,  a  lot  of  models  exist,  based  on  several  different 
hypothetical  scattering  mechanisms.  They  include  models  based  on  scattering  at 
the  sediment  interface,  like  those  of  Patterson1  and  Clay  and  Medwin,2  as  well 
as  the  volume  scattering  models  of  Nolle,3  Stockhausen,4  and  Ivakin  and 
Lysanov.5  Jackson  et  al.6-8  developed  a  two-component  model  that  treats 
interface  roughness  scattering  and  volume  scattering  separately.  Applying  this 
model  to  several  sites  suggests  that  both  components  are  important,  each 
dominating  the  total  backscatter  prediction  in  half  of  six  cases. 

Jackson’s  model  characterizes  the  seabed  with  six  parameters,  two  of 
which  describe  the  interface  roughness  and  three  of  which  describe  acoustic 
propagation  through  the  sediment,  modeled  as  an  acoustic  fluid.  The  remaining 
parameter  is  an  empirically  determined  quantity  which  specifies  the  level  of 
volume  scattering  within  the  sediment.  No  attempt  is  made  to  specify  the  actual 
volume  scattering  mechanisms  involved. 

In  1990,  Chotiros9  began  a  comprehensive  compilation  of  all  available 
shallow  grazing  angle  backscatter  data  published  in  the  literature.  The  results 
appeared  to  show  some  statistically  significant  trends  that  suggest  specific 
volume  scattering  mechanisms,  that  might  be  modeled  physically.  Two 
hypothetical  scattering  mechanisms,  involving  scattering  from  sediment  grains 
and  from  trapped  bubbles,  appear  to  fit  the  general  trends  in  the  data.  At  present 
the  data  set  is  expanding  and  the  trends  remain. 

In  1992  the  results  of  three  acoustic  penetration  experiments  at  sea  and  in 
the  laboratory  became  available.  They  suggested  that  sandy  sediments  might 
best  be  modeled  structurally  via  the  Biot  poroelastic  theory,  particularly  for  high 
frequencies  and  shallow  grazing  angles. 
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In  1993,  under  Contract  N00039-91-C-0082,  (entitled  "Bottom  Backscatter 
from  Trapped  Bubbles"),  development  of  a  new  backscatter  model  was  begun 
with  the  intention  of  taking  advantage  of  the  most  recent  experimental 
information.  It  was  continued  in  1994  under  the  same  contract.  The 
development  is  described  in  two  technical  reports  (ARL-TR-93-1510  and 
ARL-TR-94-2111).  The  model,  called  BOGGART,  is  intended  for  minehunting 
applications  and  is  therefore  optimized  for  shallow  grazing  angle,  high  frequency 
behavior.  Like  Jackson’s  model,  BOGGART  includes  volume  and  interface 
scattering  components.  It  is  different  in  its  modeling  of  the  sediment’s  structure 
and  in  the  scattering  mechanisms. 

BOGGART  is  unique  among  current  scattering  models  in  two  important 
ways.  First,  the  sediment  is  modeled  structurally  via  the  Biot  theory.  Second,  the 
volume  scattering  is  modeled  physically  rather  than  empirically.  BOGGART 
consists  of  three  modular  components  that  predict  partial  scattering  strengths 
due  to  sediment  grains,  gas  bubbles,  and  interface  roughness.  Hence  the  name 
BOGGART,  which  is  the  acronym  for  bottom  grain  gas  and  roughness  technique. 

This  report  describes  refinements  that  have  been  incorporated  in 
BOGGART  since  ARL-TR-49-21.  Both  of  these  affect  the  trapped  gas  bubble 
scattering  component  of  the  model.  They  are  (1)  a  modification  to  allow 
computation  of  bistatic  scattering  strengths,  and  (2)  depth  dependent  gas 
fractions.  Section  2  is  a  review  of  the  expression  for  the  backscattering  strength 
of  a  sediment,  due  to  trapped  gas  bubbles  in  resonance.  Section  3  describes  the 
calculation  of  bistatic  scattering  strengths.  Section  4  describes  modeling  of 
sediments  with  depth  dependent  gas  fractions. 
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2.  BACKSCATTERING  STRENGTH  OF  A  GASSY  SEDIMENT 


Figure  2.1  illustrates  a  method  of  determining  the  acoustic  pressure 
returned  to  an  acoustic  source  in  the  water  from  an  element  of  sediment  volume 
dxdydz  below  the  sediment  interface.  It  is  based  on  the  principle  of  acoustic 
reciprocity,  which  states  that,  in  a  linear  medium,  the  position  of  an  identical 
source  and  receiver  can  be  swapped  with  no  change  in  the  ratio  of  receive  to 
transmit  signals.  The  swapping  of  positions  can  be  interpreted  to  represent  the 
backscatter  case  where  a  scatterer  acts  as  an  acoustic  source  and  the  original 
projector  as  a  receiver. 

The  acoustic  source  and  scattering  elements  are  both  assumed  small  in 
comparison  with  the  acoustic  wavelength.  Both  are  surrounded  with  equal  virtual 
spheres  of  radius  much  larger  than  the  wavelength.  The  backscattering  of  sound 
from  the  volume  element  can  be  interpreted  as  follows.  First,  the  virtual  sphere 
surrounding  the  source  oscillates  with  surface  velocity  vo.  This  induces  an 
acoustic  pressure  pi  at  the  scattering  element.  The  scatterer  responds  with  its 
own  sphere's  surface  velocity  vi,  which  induces  an  acoustic  pressure  P2  back  at 
the  source  in  the  water  column.  These  pressures  and  surface  velocities  can  be 
related  by  acoustic  reciprocity: 


Pl_  P2 
V0"  V! 


(2.1) 


The  surface  velocities  vo  and  v-j  can  be  expressed  in  terms  of  local  acoustic 
pressures  and  impedances: 


v0  = 


Po 

*0 


(2.2) 


Vi 


epi 


(2.3) 
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Source 


(a)  Acoustic  propagation  of  incident  sound 


(b)  Acoustic  propagation  of  backscattered  sound 


Figure  2.1 

Determining  backscattered  pressure  from  an  element  of 
sediment  volume  by  acoustic  reciprocity. 
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where  vo,  Po.  and  Zo  are  the  surface  velocity,  acoustic  pressure,  and  acoustic 
impedance  at  the  surface  of  the  virtual  sphere  surrounding  the  source,  vi,  pi, 
and  Zi  are  the  surface  velocity,  acoustic  pressure,  and  impedance  at  the  sphere 
surrounding  the  scattered  e  is  a  transfer  function  from  incident  pressure  upon  the 
scatterer  to  scattered  pressure  at  its  virtual  sphere.  Its  average  square 
magnitude  is  expressible  in  terms  of  the  scattering  cross  section  of  the  volume 
element: 

(|e|!)=a^f  ,  (2.4) 

'  '  4jcr| 

where  rs  is  the  chosen  radius  for  the  virtual  spheres  and  ot>v  is  an  average 
scattering  cross  section  per  unit  volume  of  sediment.  The  average  is  taken  over 
the  entire  ensemble  of  possible  distributions  of  scatterers  within  the  volume 
element.  If  multiple  scattering  is  neglected,  it  can  be  expressed  as 

obv  =  Jo  N(rb)a(rb)drb  ,  (2.5) 


where  N(rt>)  is  the  bubble  size  density  function  and  rb  is  the  bubble  radius.  The 
bubble  size  density  function  is  estimated  from  the  grain  size  density  function,  as 
described  in  Ref.  1 0. 


Combination  of  Eqs.  (2.1)  -  (2.4)  yields  an  expression  for  the  pressure 
returned  to  the  projector  from  the  element  dxdydz  of  sediment  volume: 


'in  l2\ 

f|p1 

> 

-G 

CM 

Z„„ 

Hr 

J  Po 

4itrf 

Z7P’ 

2 

dxdydz 


(2.6) 


Recent  experiments12-13  suggest  that  the  best  way  to  model  the  acoustics  of 
sandy  sediments  is  via  the  Biot  theory  for  fluid  filled  porous  media.14-16  Biot's 
theory  predicts  two  compressional  waves  to  exist  in  the  pore  fluid  surrounding  the 
bubbles.  These  waves  both  contribute  to  the  total  backscattered  pressure: 
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where  Z-jf,  pit,  Zis,  and  pis  are  partial  acoustic  impedances  and  refracted 
pressures,  obtainable  via  the  Biot  theory.  Obvf  and  abvs  are  scattering  cross 
sections  that  depend  on  whether  a  Biot  fast  wave  or  slow  wave  is  being 
scattered.  These  quantities  are  described  in  more  detail  in  Ref.  1 1 . 

The  effective  scattering  Pbs  from  a  sediment  interface  can  be  regarded  as 
a  sum  of  contributions  from  all  volume  elements  below  the  interface.  In  integral 
form, 


The  backscattering  strength  of  the  interface  is  defined  as 

(|p.l2) 

BS  =  1 0  log  y - 4-  ,  (2.9) 

I P  inc  I 

where  pmc  is  the  acoustic  pressure  of  a  plane  wave  incident  upon  a  unit  area  of 
the  sediment  interface,  and  ps  is  the  pressure  of  a  scattered  wave  at  a  distance 
of  1  m  from  the  interface  element,  neglecting  absorption  and  scattering  within  the 
water  column.  The  pressures  pmc  and  ps  can  be  related  to  pressures 
experienced  at  the  projector,  by  accounting  for  spherical  spreading  and 
attenuation  in  the  water  column, 

|p8|  =  |Pbs|  r^er8<x  •  (2-10) 

|Plnc|  =  |  Po  |  T^e-rs0t  •  (2-11) 
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where  a  is  the  acoustic  absorption  in  the  water  column,  rim  is  the  unit  radius  and 
r  is  the  distance  from  the  interface  element  to  the  projector.  Substitution  of 
Eqs.  (2.8),  (2.10),  and  (2.11)  into  Eq.  (2.9)  yields  an  expression  for  the 
backscattering  strength: 


f 


\ 


BS  =  10  log 


1 

f*  M 

In  I2 

f 

^bvf 

2° 

2  _ 

■  ®bvs 

|Z0_  1 

\ 

2 

p2 

V'lmJ 

• 

|Pl| 

0 

4jt 

z7,p" 

+  4ic 

N 
«  ' 

TJ 

CO 

/ 

|Pinc| 
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•  (2.12) 


7 


This  page  intentionally  left  blank. 


8 


3.  BISTATIC  SCATTERING 


Figure  3.1  is  an  illustration  of  the  bistatic  scattering  problem.  A  projector 
at  A  insonifies  an  element  dxdydz  of  sediment  volume,  which  scatters  sound.  A 
receiver  exists  in  the  water  column  at  C.  The  principle  of  acoustic  reciprocity  is 
invoked  to  determine  a  relationship  between  the  sound  scattered  at  B  and  that 
measured  at  C. 


The  procedure  is  as  follows.  First,  the  scatterer  at  B  and  the  receiver  at  C 
are  surrounded  by  virtual  spheres  with  radius  rs  »A,.  The  pressure  at  B  that 
would  be  induced  by  a  unit  velocity  at  the  surface  of  the  virtual  sphere  at  C  can 
be  computed  by  way  of  the  Biot  theory.  By  reciprocity,  the  pressure  at  C  due  to  a 
unit  velocity  at  the  surface  of  B  is  the  same.  This  pressure  is  then  multiplied  by 
the  scattered  surface  velocity  at  B  that  is  produced  when  B  is  insonified  by  the 
source  at  A.  The  square  magnitude  of  pressure  at  the  receiver  due  to  bistatic 
scattering  from  a  volume  element  dxdydz  is  given  by 


(H’W. 


Z°n 

Z^P« 


_£o_D 

7  PlS 


dxdydz 


where  pb  is  the  total  pore  pressure  at  the  scattering  element  that  would  be 
induced  by  a  pressure  pc  at  the  surface  of  the  virtual  sphere  of  radius  rs 
surrounding  the  receiver  at  C.  Just  as  is  done  in  Eq.  (2.8)  for  the  backscatter 
case,  the  square  pressure  from  an  interface  element  dxdy  is  defined  as  the  sum 
of  contributions  from  all  volume  elements  below  dxdy, 
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Receiver 


Projector 

A 


Figure  3.1 


Calculation  of  bistatic  scattering  strength  by  reciprocity. 
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The  interface  bistatic  scattering  strength  can  be  defined  in  terms  of  the  scattered 
and  incident  pressures  at  the  interface: 


SS  =  1 0  log 


(3.3) 


where  pjnc  is  the  incident  acoustic  pressure  upon  the  interface  and  ps  is  the 
scattered  pressure  at  unit  distance  from  the  interface  element  dxdy.  ps  can  then 
be  related  to  the  pressure  pfS  experienced  at  C  by  accounting  for  absorption  and 
spherical  spreading: 


Ps  |  =  |  Pfs 


erccc 


(3.4) 


where  rc  is  the  distance  from  surface  element  dxdy  to  the  receiver.  The  bistatic 
scattering  strength  of  the  interface  is  given  by  substitution  of  Eqs.  (3.2)  and  (3.4) 
into  Eq.  (3.3): 


Figure  3.2  is  a  plot  of  the  bistatic  scattering  strength  versus  incident  and 
exit  grazing  angle.  The  modeled  sediment  is  described  in  Table  3.1.  Reflection 
from  the  water-sediment  interface  is  not  included  in  this  plot.  There  is  no 
azimuthal  dependence  since  the  bubbles  are  small  in  comparison  with  the 
acoustic  wavelength  and  are  excited  as  monopoles. 
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Scattering  strength  -  dB 


Figure  3.2 

Bistatic  scattering  strength  for  sediment  described  in  Table  3.1, 
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Table  3.1 

Bistatic  scattering  mode!  sit© 


Porosit 
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Mean  Grain  Diameter 
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Pore  Size  Parameter 
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Virtual  Mass  Parameter 


Grain  Bulk  Modulus 


Frame  Shear  Modulus 


Shear  Log  Decrement 


Frame  Bulk  Modulus 


Bulk  Log  Decrement 


Gas  Bulk  Modulus 


Gas  Densit 


Gas  Heat  Conductivit 


Gas  Specific  Heat  (const  press' 


Gas  Specific  Heat  Ratio,  Cp/Cv 


Bubble  Surface  Tension 


Bubble/Pore  Radius  Ratio 


Gas  Content 


*  Parameters  apply  to  a  sandy  site  near  Panama  City,  Florida20 


fractions  have  been  indicated  elsewhere,  including  a  laboratory  sediment  tank  at 
ARL:UT,  where  an  absence  of  gas  in  the  first  few  centimeters  of  sediment  has 
been  observed.  Hence  the  motivation  for  a  revised  version  of  Eq.  (2.12)  that 
allows  for  depth  dependent  gas  fractions. 

Two  of  the  components  of  Eq.  (2.12)  that  will  be  affected  by  a  varying  gas 
fraction  are  the  volume  scattering  cross  sections  obvf  and  abvs-  Each  is 
proportional  to  the  density  of  scatterers  and  therefore  to  the  gas  fraction.  In  the 
absence  of  multiple  scattering,  they  are  expressed  as  integrals  over  all  possible 
bubble  radii: 


N(rb)  of(rb)drb 


N(rb)  os(rb)drb 


where  rb  is  the  bubble  radius  and  N(rb)  is  the  density  function  of  bubble  sizes, 
which  is  defined  as  the  number  of  bubbles  per  unit  volume  per  unit  radius 
increment,  at  and  as  are  single  bubble  scattering  cross  sections  for  fast  and  slow 
waves,  respectively.  They  are  given  by  Wildt:18 


of(rb)  = 


4  -1  +5, 2 
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trapped  gas,  which  scatters  and  absorbs  sound.  The  fast  and  slow  wave 
pressures  Pbf(z)  and  Pbs(z)  in  a  gassy  sediment  can  be  expressed  as 


Pbf(Z)  =  Pbfo(Z)  m 

(4.9) 

Pbs(Z)  ”  Pfoso(Z)  P  S(Z)  > 

(4.10) 

where  Pbfo(z)  and  Pbso(z)  are  the  fast  and  slow  wave  acoustic  pressures  that 
would  exist  if  no  gas  were  present.  (3f(z)  and  ps(z)  are  multipliers  that  adjust  for 
scattering  and  absorption  of  sound  by  trapped  gas.  They  can  be  expressed  in 
terms  of  the  sediment’s  depth  dependent  extinction  cross  section, 


P  f(z)  =  exp 


fZ^fev(Zi) 

Jo  sinO 


(4.11) 


P  s(z)  =  exp 


*z<WZi) 

o  sin  O 


dzi 


(4.12) 


where  is  the  grazing  angle.  The  fast  and  slow  wave  sediment  extinction  cross 
sections  otev(z)  and  asev(z)  can  be  expressed  in  terms  of  the  gas  fraction  and  the 
bubble  size  density  function, 


<WZi)  =  C(Zi)  J~F(rb)  afe(rb)  drb  ,  (4.13) 

°sev(Zi)  -  C(Z-j)  F(^b)  ^se(^b)  ^ b  >  (4. 1 4) 

where  £(zi)  is  the  depth  dependent  gas  fraction  and  F(rb)  is  the  bubble  size 
density  function.  Ofe(rb)  and  ose(rb)  are  the  fast  and  slow  wave  extinction  cross 
sections  of  a  single  bubble  of  radius  rb,  given  by  Wildt.18 

°fe(rb)  =  Sf(rb)  ^  .  (4.15) 
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(4.16) 


where  Gf(rb)  and  os(rb)  are  the  fast  and  slow  wave  scattering  cross  sections,  8f 
and  Ss  are  the  corresponding  damping  constants,  and  k  is  the  acoustic 
wavenumber. 

Figure  4.1  illustrates  the  effect  of  depth  dependent  gas  fractions  on  the 
backscattering  strength.  It  is  a  comparison  of  predicted  backscattering  strength 
versus  grazing  angle  over  a  sandy  sediment  bottom  with  different  gas  fraction 
profiles.  The  sediment  geoacoustic  parameters  are  those  listed  in 
Table  4.1.  The  gas  fraction  profiles  include  (1)  an  infinite  homogeneous  half 
space  with  gas  fraction  of  1 0'5,  (2)  a  7  cm  gasless  layer  over  an  infinite  half 
space  with  gas  fraction  of  1 0-5,  and  (3)  a  20  cm  gasless  layer  over  an  infinite  half 
space  with  gas  fraction  of  1 0'5. 

In  Fig.  4.1,  the  predicted  backscattering  strengths  have  been  shifted 
vertically  to  line  up  the  initial  portions  of  the  scattering  strength  curves.  It  is 
apparent  from  this  figure  that  in  addition  to  decreasing  the  backscattering 
strength,  the  addition  of  a  gasless  layer  has  the  effect  of  causing  the 
backscattering  strength  to  increase  more  sharply  with  grazing  angle.  The  dashed 
curve  in  Fig.  4.1  is  a  plot  of  Lambert’s  rule,  defined  as 

BS  =  20  log(sinQ)  +  r|  ,  (4.17) 

where  0  is  the  grazing  angle  and  rj  is  a  constant,  which  can  be  determined  by  a 
fit  to  experimental  data.  Lambert’s  rule  has  been  widely  observed  to  fit  closely  to 
experimental  data,20  but  has  no  rigorous  theoretical  grounding.  In  Fig.  4.1, 
Lambert’s  rule  appears  to  follow  the  7  cm  layer  data  most  closely.  This  suggests 
that,  in  the  case  of  scattering  by  sediment  bubbles,  Lambert’s  rule  scattering  may 
be  a  manifestation  of  depth  dependent  scattering  properties. 
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Figure  4.1 

Baekscattering  strength 

versus  grazing  angle  for  different  gas  fraction  profiles. 
Sediment  parameters  are  given  in  Table  4.1. 
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5.  CONCLUSIONS 

The  ARL:UT  bottom  backscatter  model,  BOGGART,  has  been  improved  in 
two  ways.  Both  of  these  improvements  affect  the  trapped  bubble  scattering 
component  of  the  model. 

The  model  has  been  generalized  for  bistatic  scattering  from  trapped  gas 
bubbles.  The  technique  is  derived  from  the  principle  of  acoustic  reciprocity, 
which  depends  on  an  assumption  that  scatterers  are  small  in  comparison  with  the 
acoustic  wavelength.  The  resulting  bistatic  model  has  no  azimuthal  dependence. 
It  is  possible  that  a  more  complete  bistatic  scattering  model,  including  sediment 
grain  and  roughness  components,  will  have  some  azimuthal  dependence. 

Depth  dependent  gas  fractions  have  been  incorporated  in  the  model. 
Comparisons  of  different  gas  fraction  profiles  suggest  that  the  grazing  angle 
dependence  of  the  backscattering  strength  is  strongly  affected.  At  present  the 
gas  fraction  is  the  only  parameter  in  the  model  that  is  allowed  to  vary  with  depth. 
A  more  complete  model  would  allow  for  depth  dependence  of  other  important 
components,  such  as  the  bubble  size  density  function  and  the  Biot  parameters. 

The  current  improvements  will  give  BOGGART  greater  applicability  to 
actual  sediments,  in  which  horizontal  layering  is  probably  common.  The 
modification  for  bistatic  scattering  will  help  a  user  make  sonar  performance 
predictions  in  a  multipath  environment.  Though  BOGGART  is  complete  and 
functional,  it  is  modular  and  designed  with  future  upgrades  in  mind. 
Improvements  will  continue  as  more  experimental  information  about  the  seafloor 
becomes  available. 
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